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Abstract 

A measurement is presented of elastic deeply virtual Compton scattering 7*p — > 7p 
made using e+p collision data corresponding to a luminosity of 46.5 pb^^, taken with 
the HI detector at HERA. The cross section is measured as a function of the photon vir- 
tuality, Q^, the invariant mass of the 7*p system, W, and for the first time, differen- 
tially in the squared momentum transfer at the proton vertex, t, in the kinematic range 
2 < < SOGeV^, 30 < 1^ < 140GeVand[t[ < 1 GeV^. QCD based calcula- 
tions at next-to-leading order using generalized parton distributions can describe the data, 
as can colour dipole model predictions. 
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1 Introduction 



Measurements of the deep-inelastic scattering (DIS) of leptons and nucleons allow the ex- 
traction of Parton Distribution Functions (PDFs) which describe the longitudinal momentum 
carried by the quarks, anti-quarks and gluons that make up the fast-moving nucleons. While 
these PDFs provide crucial input to perturbative Quantum Chromodynamic (QCD) calculations 
of processes involving hadrons, they do not provide a complete picture of the partonic struc- 
ture of nucleons. In particular, PDFs contain neither information on the correlations between 
partons nor on their transverse motion. This missing information can be provided by measure- 
ments of processes in which the nucleon remains intact, such as the exclusive production of light 
meson states in lepton-nucleon collisions, and is encoded in Generalised Parton Distributions 
(GPDs) [1-4]. 

The simplest process sensitive to GPDs is deeply virtual Compton scattering (DVCS) (fig- 
ure[Tt), which is the diffractive scattering of a virtual photon off a proton [5-10], 7*p 7p. In 
the present analysis DVCS is accessed through the reaction: 

e~^p — > e'^jp. (1) 

This process is of particular interest as it has both a clear experimental signature and is 
calculable in perturbative QCD: it does not suffer from the uncertainties caused by the lack of 
understanding of the meson wave function that plague exclusive vector meson electroproduc- 
tion. 




Figure 1: Diagrams illustrating the DVCS (a) and the Bethe-Heitler (b and c) processes. 

The reaction studied receives contributions from both the DVCS process, whose origin lies 
in the strong interaction, and the purely electromagnetic Bethe-Heitler (BH) process (figures [1J5 
and[lj;), where the photon is emitted from the positron. The BH cross section can be precisely 
calculated in QED using elastic proton form factors. Here, the DVCS cross section is obtained 
by subtracting the BH contribution from the total cross section, which is possible since the 
interference contribution vanishes [10], as this measurement is integrated over azimuthal angles. 

The first measurements of the DVCS cross section at high energy were obtained by HI [11] 
and ZEUS [12] and the helicity asymmetry in DVCS has been measured at lower energy with 
polarised lepton beams by HERMES [13] and CLAS [14]. 
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In this paper, a measurement of the DVCS cross section is presented, based on data collected 
with the HI detector at HERA in the years 1996 to 2000. These data correspond to a luminosity 
of 46.5 pb^^, a factor of 4 larger than the luminosity used in the previous HI publication [11], 
which is based only on 1997 data. The cross section is presented as a function of the photon 
virtuality, , the invariant mass of the 7*p system, W , and the squared momentum transfer at 
the proton vertex, t. 

2 Generalized Parton Distributions and Theoretical Predic- 
tions 

The leading order diagram for DVCS in positron proton scattering is shown in figure and 
a diagram that contributes at next-to-leading order in figure |2h- The transition from a virtual 
photon to a real photon forces the fractional momenta of the two partons involved to be dif- 
ferent ("skewed"). Hence, DVCS is sensitive to the correlations between partons in the proton 
which are encoded in the GPDs. In the presence of a hard scale, here Q^, the DVCS scattering 
amplitude factorises [3,6,7] into a hard part, calculable order by order in perturbative QCD, and 
the GPDs which contain the non-perturbative effects due to the structure of the proton. 



Figure 2: Examples of diagrams for the DVCS process a) at leading order, b) at next-to-leading 
order 



2.1 Generalized Parton Distributions 

The GPDs generalize and interpolate between the PDFs and elastic form factors. The PDFs 
contain information on the longitudinal momenta of the partons while form factors contain 
information on their transverse momenta, often in the form of sum rules related to charges, local 
currents and the energy-momentum tensor of QCD. GPDs have simple physical significance in 
light-cone coordinates (or the infinite momentum frame), where they represent the interference 
of two different wave functions, one of a parton having a momentum fraction x + ^ and the 
other of a parton with a momentum fraction x — ^, as is illustrated in figure |2| Besides the 
longitudinal momentum fraction variables ^ (called skewedness) and x, GPDs depend on t, the 



e 




a) 



b) 
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square of the four-momentum exchanged at the hadron vertex. GPDs are defined at a starting 
scale //^ and their evolution is generated by perturbative QCD. 

There are two different types of GPDs (for a quark g or a gluon g) in the unpolarised case: 

H'^'^{x, ^, t) and E'''^{x, ^, t). While the E'^'^ distributions have no equivalent in the ordinary 
PDF approach, the H'^'^ reduce to the usual PDFs in the forward limit = 0, t = 0), 
i.e. H'^{x, 0, 0) = q{x) and H^{x, 0, 0) = xg{x), where q{x) and g{x) are the ordinary parton 
distributions. The variable x is defined in the range [—1,1], with negative values corresponding 
to anti-quark distributions: H'^{—x, 0, 0) = —q{x). The gluon GPD is symmetric in x in the 
forward limit: H^{—x, 0, 0) = H^(x, 0, 0). The skewedness variable ^ is related to the well 
known Bjorken-,x variable, XBj,i>y ^ = axBj/{2 — axsj), where a = 1 + q'^/Q^ and q' 
denotes the four-momentum of the outgoing photon^ The first moments of the GPDs in x are 
given by form factors [4]. 

Two different kinematic regions exist for GPDs with respect to the variables x and ^. The 
DGLAP region, where \x\ > { [15-18], corresponds to the emission and re-absorption of a 
quark, anti-quark or a gluon. The ERBL [19,20] region, where \x\ < ^, corresponds to meson 
or gluon pair exchange. Each region has its own evolution equations. 

The recent strong interest in GPDs was stimulated by the information they contain on the 
spin structure of the nucleon. In particular, GPDs are so far the only known means of probing the 
orbital motion of partons in the nucleon through Ji's Sum Rule [21], which relates unpolarised 
GPDs to the total angular momentum of the proton. DVCS measurements at HERA can provide 
constraints on this sum rule through their sensitivity to the GPDs. 

2.2 Theoretical Predictions 

The measurements presented here are compared with NLO QCD calculations and predictions 
made using colour dipole approaches. In NLO QCD, the DVCS cross section has been calcu- 
lated [22,23] using two different GPD parameterisations [24]. The t dependence of the GPDs 
is taken to be e'^^^l The MRST2001 [25] and CTEQ6 [26] parameterisations of the PDFs are 
used in the DGLAP region > 0. Thus H, which provides the main contribution to DVCS 
at small xbj, is given at the starting scale /i by H'^{x,^,t) = q{x) e"*'*' for the quarks and 
H^(x, ^,t) = X g{x) e"''!*' for the gluons^. Both the skewing and the Q'^ dependence are gener- 
ated dynamically. In the ERBL region {\x\ < 0, these parameterisations have to be modified, 
ensuring a smooth continuation to the DGLAP region (for details see [24]). These GPD models 
are found to describe both the shape of the previous HI DVCS cross section measurements [11] 
and the single spin asymmetry measured by HERMES [13]. 

The DVCS cross section has also been calculated in the colour dipole approach, which 
is successful in describing both inclusive and diffractive scattering in the DIS regime at high 
energy. These predictions are based on a factorisation of the DVCS amplitude into the wave 
function for the photon to fluctuate into a qq pair, the cross section for this pair to interact with 

'For the DVCS process, the outgoing photon is real {q'"^ = 0) and ^ reduces to XBj/{'i — XBj)- The forward 
limit corresponds to the case of inclusive DIS, where g'^ = —Q"^ and thus ^ = 0. 

different ansatz for GPDs has been used in [27] in a LO calculation of the DVCS cross section. 
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the proton and the outgoing photon wave function. If s-channel helicity is conserved in DVCS, 

the virtual photon must be transversely polarised. As the wave function of the transversely 
polarised 7* can select large dipole sizes, whose interactions are predominantly soft, DVCS 
constitutes a good probe of the transition between the perturbative and non-perturbative regimes 
of QCD. The various calculations differ in the way the dipole cross section is parameterised. 
Donnachie and Dosch [28] use soft and hard pomeron exchange depending on the size of the 
dipole. All parameters are determined fmmpp and total cross section measurements. Favart 
and Machado [29] apply the saturation model of Golec-Biemat et al. [30] to the DVCS process 
and use DGLAP evolution [31], following the approach of Bartels, Golec-Biernat and Kowalski 
(BGBK) [32]. In both cases an exponential ^-dependence, e~^l*l, is assumed. 

3 Experimental Procedure 

3.1 HI Detector 

A detailed description of the HI detector can be found in [33]. Here only the detector com- 
ponents relevant for the present analysis are described. The SpaCal [34], a lead scintillating 
fibre calorimeter, covers the backward^ region of the HI detector (153° < 9 < 177.5°). Its en- 
ergy resolution for electromagnetic showers is a{E)/E ~ 7.1%/\/E/GcV © 1%. The liquid 
argon (LAr) calorimeter (4° < 6^ < 154°) is situated inside a solenoidal magnet. The energy 
resolution for electromagnetic showers is a {E)/E ~ 11%/v/^/GeV as obtained from test 
beam measurements [35]. The backward drift chamber (BDC), placed in front of the SpaCal, 
measures track segments for charged particles entering the SpaCal from the interaction region. 
These are used to identify the scattered positron and to determine its position with a resolution 
of 0.5 mm in the radial and 2.5 mm in the azimuthal direction. The main component of the 
central tracking detector is the central jet chamber (CJC) which consists of two 2 m long coax- 
ial cylindrical drift chambers, with wires parallel to the beam direction. The measurement of 
charged particle transverse momenta is performed in a magnetic field of 1.15 T, uniform over 
the full tracker volume. The forward components of the detector, used here to tag hadronic ac- 
tivity at large pseudo-rapidity (5 £ 7), are the forward muon detector (FMD) and the proton 
remnant tagger (PRT). The FMD, designed to identify muons emitted in the forward direction, 
contains six planes of drift cells. It is used here to detect the particles produced when a proton 
dissociates and secondary interactions occur in the beampipe and adjacent material. Secondary 
particles, or the scattered proton, can also be detected by the PRT, which is located at 24 m 
from the interaction point and consists of layers of scintillator surrounding the beam pipe. The 
luminosity is determined from the rate of BH events measured in a luminosity monitor. 

3.2 Kinematics 

For DVCS, the final state photon does not originate from the positron and therefore the ratio of 
the DVCS to the BH cross sections is expected to increase when the photon is scattered in the 

^Hl uses a right-handed coordinate system with z axis along the beam direction, the +z or "forward" direction 
being that of the outgoing proton beam. The polar angle is defined with respect to the z axis and the pseudo- 
rapidity is given by 77 = — In tan 0/2. 
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forward direction. The analysis sample is thus selected by requiring a photon candidate in the 
LAr calorimeter and a positron candidate in the SpaCal calorimeter. 

The reconstruction of the kinematic variables Q^, xbj and W relies on the polar angle 
measurements of the final state positron, Of., and photon, 9^: 



g2 _ sin (1 + cosge) ^2) 

° sin 9-y + sin 9e — sin (6'e + 9^) ' 
^ _ Eq sin 9y + sin 9e + sin {9e + 9^) 
Ep sin + sin 9^ — sin (6'e +9^) 

W = (4) 

XBj 

where Eq and i?p are the positron and proton beam energies, respectively. For the majority of the 
events, the scattered positron trajectory is not measured in the CJC and the event vertex cannot 
be determined. The polar angles of the positron and photon are then determined assuming that 
they come from the nominal event vertex. The square of the four-momentum transfer to the 
proton, t, is very well approximated by the square of the vector sum of the transverse momenta 
of the final state photon, , and of the scattered positron, pt^ : 

(5) 



3.3 Monte Carlo Simulation 

Monte Carlo (MC) simulations are used to estimate the corrections that must be applied to 
the data due to the finite acceptance and resolution of the detector. Elastic DVCS events in 
ep collisions are generated using the Monte Carlo generator MILOU [36], which is based on 
a NLO QCD cross section calculation [22, 23, 37] (see section IT^ . and using a slope in t of 
6 = 6 GeV^. Higher order photon radiation from the incoming positron is implemented in the 
coUinear approximation. DVCS events in which the proton dissociates into a baryonic system 
Y are also simulated with the program MILOU using a t slope of bpdiss = 1-5 GeV^^ [38]. The 
Monte Carlo generator COMPTON 2.1 [39,40] is used to simulate both elastic and inelastic 
BH events. Hadronisation processes in inelastic BH events are simulated using the SOPHIA 
model [41]. Diffractive u and (p meson events are generated with the DIFFVM Monte Carlo 
program [42]. The events generated using all these programs are passed through a detailed 
simulation of the HI detector and are subject to the same reconstruction and analysis chain as 
the data. 



3.4 Event Selection 

The data were obtained with the HI detector when the HERA collider was operated with 
820 GeV (1996-1997) and 920 GeV (1999-2000) protons and 27.6 GeV positron beams. The 
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data sample corresponds to an integrated luminosity of 46.5 pb~^, 11.5 pb^^ of which were 
accumulated in 1996-1997 and 35 pb^^ in 1999-2000. The event trigger used is based on the 
detection of an energy deposition greater than 6 GeV in the electromagnetic section of the 
SpaCal calorimeter. Due to the different trigger settings, selected events in the 1996-1997 pe- 
riod are in the kinematic range > 2 GeV^ while those in the 1999-2000 period are in the 
range g2 > 4 QeV^. 

The DVCS event selection requires that the following criteria be fulfilled. The scattered 
positron must be detected in the SpaCal, have an energy larger than 15 GeV and be validated 
by a track segment in the BDC. The photon must be measured in the LAr calorimeter with a 
transverse momentum > 1 GeV (1996-1997) orp^ > 1.5 GeV (1999-2000) and a polar angle 
between 25° and 145°. The scattered proton escapes undetected through the beam pipe. Events 
with more than one central track are rejected while events with one central track are only kept if 
that track is associated with the scattered positron. In order to reject inelastic and proton disso- 
ciation events, no further energy deposition in the LAr calorimeter with energy above 0.5 GeV 
is allowed and no activity above the noise level is allowed in the PRT and FMD. The influence 
of QED radiative corrections is reduced by the requirement that the longitudinal momentum 
balance — Pz) > 45 GeV. Here, E denotes the energy and Pz the momentum along the 
beam axis of the final state particles and the sum runs over all such particles. To enhance the 
DVCS signal with respect to the BH contribution and to ensure a large acceptance, the kinematic 
domain is explicitly restricted to < 80GeVMt| < 1 GeV^ and 30 < W < 140 GeV. 

The selected sample contains 1243 events and is dominated by the DVCS contribution, but 
also contains contributions from the elastic BH process and from the (inelastic) BH and DVCS 
processes with proton dissociation, e^p e^jY, where the baryonic system Y of mass My is 
not detected in the forward detectors. 

As in previous HI DVCS analyses [11,43], a control sample of BH events is also selected. 
Here, it is required that the positron be detected in the LAr and the photon in the SpaCal. It has 
been verified that the COMPTON MC correctly describes the normalisation and the shapes of 
the distributions of the kinematic variables for these events within an uncertainty of 5%. Using 
events with a signal in the forward detectors, and subtracting the inelastic BH contribution, 
obtained from the COMPTON MC, the contribution of proton dissociation to the DVCS event 
sample is estimated to be 16 ± 8 % for the 1996-1997 data (lower Q^) and 10 ± 5 % for the 
1999-2000 data. The other backgrounds considered are diffractive cu and (f) production, with 
decay modes to final states including photons. The main backgrounds originate from the decays 
to 7r°7 and ^l^s followed by the decay Kg vr^vr^. The contribution of these 
processes to the DVCS sample is estimated to be below 3.5% for the data taken in 1996-1997 
and below 1% for that taken in 1999-2000. 

In figure 13 the data are compared with the sum of the MC expectations. The BH contribu- 
tions and the u and (p backgrounds are normalised to the luminosity. The DVCS contribution 
is normalised such that the sum of the DVCS, BH and diffractive vector meson contributions is 
equal to the total number of events in the data. The distributions of the energy and polar angle 
of the positron and the photon are shown in figures |3^-d. The coplanarity, shown in figure|3^, is 
defined to be the difference of the azimuthal angles of the electron and photon directions. It is 
related to the pf-balance of the positron-photon system. The distribution of the invariant mass 
of the positron and the photon is presented in figure |3f. The sum of the MC contributions gives 
a good description of the shapes of the data distributions. 
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Figure 3: Distributions of the energy of the scattered positron (a), the energy of the photon (b), 
the polar angle of the scattered positron (c), the polar angle of the photon (d), the coplanarity 
(e) and the positron-photon invariant mass (f). The data are compared with MC expectations 
for elastic DVCS, elastic BH, BH and DVCS with proton dissociation, and to and (p diffractive 
backgrounds. The DVCS contribution is normalised such that the sum of the DVCS, BH and 
diffractive vector meson contributions is equal to the total number of events in the data. The 
normalisation of the other contributions is described in the text. 
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3.5 Cross Section Measurement Method 



To extract the cross section, the selected data are corrected for detector efficiencies, acceptance, 
bin-to-bin migrations and for initial state radiation from the positron using the Monte Carlo 
simulation. The inelastic BH contribution is subtracted bin by bin using the COMPTON Monte 
Carlo program. The contribution of DVCS events with proton dissociation is subtracted bin by 
bin using the MILOU Monte Carlo simulation. A 5% correction is applied to correct for the loss 
of elastic DVCS events due to the requirement that there be no signal in the forward detectors. 
The background contributions from diffractive uj and production are also subtracted using the 
MC simulations. 

In the leading twist approximation, the main contribution resulting from the interference of 
the BH and DVCS processes is proportional to the cosine of the azimuthal angle of the photon'*. 
Since the present measurement is integrated over this angle, the overall contribution of the 
interference term is negligible. The elastic BH cross section can therefore be subtracted from 
the total e^p e^'jp cross section in order to obtain the contribution from DVCS processes. 
This contribution is then converted to the ^p cross section using the equivalent photon 

approximation^: 



where the transverse photon flux T is given by [44], 



2 



Here, s is the square of the ep centre-of-mass energy. 
The t dependence is factorised according to: 



da[Yp ip] da[7> 7p] 



with y = . (7) 



dt ^ 



e-''l*l . (8) 

i=0 



The cross section a['j*p ^p] is extracted from equations|3and[H]using an iterative proce- 
dure and fitting the t integrated cross section with the form: 

c^b*P^lp]{Q\y)=N-y'/'-(^-^^ , (9) 

where 5, n and h are free parameters and N is fixed by the integration of equation |6l More 
details can be found in [43]. 

The same method is used to extract a[y*p 7p] as a function of and of xsj- 



^The azimuthal angle of the photon is defined as the angle between the plane formed by the incoming and 
scattered positron and that formed by the 7* and the scattered proton. 

^After integrating over azimuthal angles only transversely polarised 7* contribute to the DVCS process. 
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3.6 Systematic Errors 



The main sources of systematic errors and their resulting uncertainty on the DVCS cross section 
measurements are: 

• the subtraction of the DVCS proton dissociation background (typically 11% in 1996- 
1997, 8% in 1999-2000 and up to 20% in the highest \t\ bin) estimated using MC simula- 
tions with fep^j^s = 1.5 ± 0.5 GeV"2 and an My dependence rfa/fiM-^ ~ (l/My)2-o±o-3 ; 

• the uncertainty on the acceptance correction factors (typically 10% and up to 25% in the 
highest \t \ bin) calculated by varying b between 4 and 7 GeV^^; 

• the uncertainty on the determination of 5 and n used for the bin centre corrections (which 
ranges between 9 and 16%); 

• the uncertainty on the BH subtraction (up to 7% for the highest W bin); 

• the uncertainties on the vertex position and the measurement of the scattered positron/photon 
angles (each contribution leading to up to 12% in the highest |t| bin); 

• the uncertainties on the positron/photon energies (each contribution leading to up to 12% 
in the highest |t| bin); 

• the noise in the CJC (typically 4%), and in the FMD (up to 2%); 

• the luminosity measurement (typically 2.5%). 

The total systematic error is found to be typically 25%. 



4 Results 



4.1 Cross Sections 

The cross sections are determined separately for the two data taking periods, which cover dif- 
ferent ranges in Q^, and are then combined. The 1996-1997 period covers the kinematic range 
2 < g2 < 20GeV2and30 < W < 120 GeV, the 1999-2000 period 4 < < 80 GeV^ 
and 30 < W < 140 GeV; in both cases |t| < 1 GeVl 

The 7*p cross section is shown differentially in t in figure |4l and given in table [T] for 
Q2 = 4 GeV2 and W = 71 GeV (using the 1996-1997 data) and = 8 GeV^ and 
W = 82 GeV (using the 1999-2000 data). The data points are fitted with the exponential form 
e"^l*l, which gives 6 = 6.66 ± 0.54 ± 0.43 GeV^^ at = 4 GeV^ where the first error is sta- 
tistical and the second systematic. At = 8 GeV^, a value of 6 = 5.82 ± 0.59 ± 0.50 GeV"^ 
is obtained. The two cross sections are averaged after correcting the 1996-1997 results to 
Q"^ = 8 GeV^ and W = 82 GeV using equation |9l (see table [B- The t slope is then 
measured to be 6 = 6.02 ± 0.35 ± 0.39 GeV^^. 
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The cross section as a function of is shown in figure |5] and given in table El for W = 
82 GeV and |t| < 1 GeV^. Fitting the Q"^ dependence with the form (l/Q^)" gives n = 
1.54 ± 0.09 ± 0.04. The Q'^ dependence of the cross section is also given for a fixed value of 
XBj = 1-8 ■ 10~^ in table|2l in the restricted range accessible for fixed xbj- 

The cross section as a function of W is shown in figure |^ and given in table |5] for = 
4 GeV^ and Q'^ = 8 GeV^; in both cases \t\ < 1 GeV^. The data are fitted using the form 
IV-^ which gives 5 = 0.69 ± 0.32 ± 0.17 at = 4 GeV^ and 5 = 0.81 ± 0.34 ± 0.22 
at = 8 GeV^. The two measurements are combined as explained above and the resulting 
cross section is given in table |3l at Q"^ = 8 GeV^. Fitting the combined sample with the form 
VT'^ gives 5 = 0.77 ± 0.23 ± 0.19. The steep rise of the cross section with is a strong 
indication of the presence of a hard scattering process, the value of 6 being comparable to that 
measured in exclusive J/^ production [45,46]. 

The extracted values of b, 5 and n are summarised in tablelH 
4.2 Discussion 

The cross section measurements from the combined data sample are shown with ZEUS mea- 
surements^ [12] and theoretical predictions as a function of Q'^ in figure |7^ and as a function 
of W in figure[8^. All predictions are made assuming an exponential dependence on \t\, using 
the measured value b = 6.02 ± 0.52 GeV^^. The error represents the total uncertainty of 
the t slope which is reflected in the band associated with each of the predicted curves. The HI 
and ZEUS measurements are seen to be consistent. The NLO QCD calculations of Freund et 
al. use two different GPDs, based on MRST 2001 and CTEQ6, for the diagonal distributions 
in the DGLAP domain. These two parameterisations show similar behaviour in and in W 
and differ mainly in the normalisation, which reflects the relative size of the quark singlet and 
gluon distributions for each set. The HI data are better described by the parameterisation based 
on CTEQ6, but it must be noted that the prediction also depends on the parameterisation of the 
ERBL region. 

As shown in figures 7b and 8b, colour dipole models also provide a reasonable description 
of the data, both in shape and in normalisation. The dependence is better described by 
the Favart-Machado prediction when DGLAP evolution of the dipole (BGBK) is included. As 
regards the W dependence, the HI data are consistent with both the Donnachie-Dosch and 
the Favart-Machado predictions, while the ZEUS measurements slightly favour the Donnachie- 
Dosch prediction. 

Introducing a dependence of the \t\ slope , b = bo{l - 0.15 log((5V2)) GeV"^ [23], 
as extracted for exclusive p meson production [47,48] (with bo such that b = 6.02 GeV~^ at 
= 8 GeV^), does not significantly change the above conclusions. 

^The ZEUS measurements, for = 89 GeV, have been rescaled to = 82 GeV and from = 9.6 GeV^ 
to = 8 GeV^ using the parameter values (5 = 0.75 and rt = 1.54 as quoted by ZEUS. 
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5 Conclusion 



The DVCS process has been studied in the kinematic region 30 < W < 140 GeV, 2 < 
Q'^ < 80 GeV^ and \t\ < 1 GeV^ using data taken with the HI detector in the years 1996 to 
2000. The 7*p cross section has been measured as a function of and as a function 

of W, and for the first time differentially in t. The dependence of the cross section on Q ^ is 
well reproduced by the shape (l/Q^)" with n = 1.54 ± 0.09 ± 0.04 at = 82 GeV. The 
W dependence can be described by a fit of the form yielding 5 = 0.77 ± 0.23 ± 0.19 
at = 8 GeV^. The fall of the cross section differential in t can be described by the form 
e-^l*l with b = 6.02 ± 0.35 ± 0.39 GeV'^ atQ^ = 8 GeV^. This first measurement of the 
t dependence of DVCS constrains the normalisation of the theoretical predictions. NLO QCD 
calculations give a good description of the normalisation as well as of the Q'^ and W dependence 
of the measured cross section using a parameterisation of the GPDs based on the CTEQ6 parton 
distribution functions. The calculations rely on ordinary (unskewed) parton distributions in the 
DGLAP region and generate the skewedness dynamically. Colour dipole model predictions also 
give a good general description of the data. This is particularly true for a saturation model in 
which the DGLAP equation is used to describe the evolution of the dipole. 
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da (7*p 7p) /dt nb/Ge 




1996-1997 


1999-2000 


All data 




g2 = 4 GeV2 


g2 = 8 GeV2 


g2 = 8 GeV2 


\t\ [GeV J 


W = 71 GeV 


W = 82 GeV 


W = 82 GeV 


0.1 
0.3 
0.5 
0.8 


29.9 ±4.1 ±7.1 
8.0 ±1.4 ±1.4 
2.13 ±0.60 ±0.69 
0.27 ±0.12 ±0.14 


13.3 ±1.9 ±3.4 
3.99 ±0.57 ±0.69 
0.90 ±0.25 ±0.30 
0.36 ±0.09 ±0.14 


12.0 ±1.2 ±2.9 
3.44 ±0.38 ±0.61 
0.84 ±0.17 ±0.29 
0.21 ±0.04 ±0.09 



Table 1: Cross sections differential in t for the two data samples and for the combined sample. 
The first errors are statistical, the second systematic. 









(t(7*p ^ 


7p) [nb] 








W 


= 82 GeV 


Xbj 


= 1.8- 


10-3 


3.0 


15.7 


±2.5 


±3.4 








5.25 


5.7 


±1.1 


±1.4 


6.74 


±0.93 


±1.02 


8.75 


3.20 


±0.49 


±0.69 


3.25 


±0.51 


±0.60 


15.5 


1.20 


±0.22 


±0.32 


1.45 


±0.30 


±0.36 


25.0 


0.70 


±0.19 


±0.19 








55.0 


0.15 


±0.05 


±0.05 









Table 2: The j*p cross section as a function of for \t\ < 1 GeV^, atW = 82 GeV 

(second column) and at xbj — 1-8 • 10-^ (third column). The first errors are statistical, the 
second systematic. 





(t(7*p — s> 7p) [nb] 


1996-1997 


1999-2000 


All data 


W [GeV] 


Q'- = 1 GeV2 


Q2 = 8 GeV2 


Q'- = 8 GeV2 


45 

70 

90 

110 

130 


6.5 ±0.8 ±1.1 
8.9 ±1.3 ±1.6 
11.1 ±2.2 ±2.7 
10.1 ±4.7 ±4.6 


2.56 ±0.36 ±0.32 
2.93 ±0.63 ±0.46 
4.45 ±0.83 ±0.82 

5.3 ±1.4 ±1.4 

6.4 ±2.5 ±2.7 


2.28 ±0.21 ±0.34 
2.91 ±0.35 ±0.51 
3.97 ±0.54 ±0.85 
4.4 ±1.0 ±1.5 
6.4 ±2.5 ±2.7 



Table 3: The 7*p 7p cross section as a function ofWfor\t\ < 1 GeV^ for the two data 
samples and for the combined sample. The first errors are statistical, the second systematic. 
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h [GeV-2] 


5 


n 


4GeV2 


6.66 ± 0.54 ± 0.43 


0.69 ±0.32 ±0.17 


1.54 ±0.09 ±0.04 


8GeV" 


5.82 ± 0.59 ± 0.50 


0.81 ± 0.34 ± 0.22 


All data, 8 GeV^ 


6.02 ± 0.35 ± 0.39 


0.77 ±0.23 ±0.19 



Table 4: Summary of the h, 5 and n values separately for the two data taking periods atQ^ — A 
GeV^ and = 8 GeV^ and for the combined sample at = 8 GeV^. The first errors are 
statistical, the second systematic. The values ofb are measured atW = 71 GeV forQ^ = 4 and 
W = 82 GeV for Q"^ = 8. The values of 6 and n are given for \t\ < 1 GeV^. The value of n is 
calculated at W ^82 GeV. 
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Figure 4: The cross section YP IP differential in t, for Q"^ = A GeV^ at W = 71 GeV 
and Q"^ = 8 GeV^ atW = 82 GeV. The inner error bars represent the statistical and the full 
error bars the quadratic sum of the statistical and systematic uncertainties. The lines represent 
the results of fits to the exponential form e^^'*', giving the values ofb shown in the insert (see 
tabJell). 




Figure 5: The 7*p — > 7p cross section as a function of for W = 82 GeV and \t\ < 1 GeV^. 
The inner error bars represent the statistical and the full error bars the quadratic sum of the 
statistical and systematic uncertainties. The curve is the result of a fit to the form 
giving the value of n shown in the figure (see table^. 
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Figure 6: The 7*p — *• cross section as a function ofW for \t\ < 1 GeV^ at = A GeV^ 
and at = 8 GeV^. The inner error bars represent the statistical and the full error bars the 
statistical and systematic uncertainties added in quadrature. The lines are the results of a fit to 
the form W^, giving the values of 5 shown in the insert (see table\^. 
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Figure 7: The 7*p 'jp cross section as a function of for W = 82 GeV and \t\ < 1 GeV^. 
The inner error bars represent the statistical and the full error bars the statistical and systematic 
uncertainties added in quadrature. The HI measurement is shown together with the results of 
ZEUS [12] and several theoretical predictions, a) Comparison with QCD predictions calculated 
at NLO by Freund et al. [24] based on MRST 2001 and CTEQ6 PDFs. b) Comparison with 
the colour dipole predictions of Donnachie and Dosch [28] and Favart and Machado with [31] 
and without [29] the DGLAP evolution of the saturating dipole (indicated as BGBK). The band 
associated with each prediction corresponds to the uncertainty on the measured t-slope. 
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Figure 8: The 7*j9 7p cross section as a function ofW for Q"^ = 8 GeV^ and \t\ < 1 GeV^. 
The inner error bars represent the statistical and the full error bars the statistical and systematic 
uncertainties added in quadrature. The measurement is shown with the results of ZEUS [12] and 
several theoretical predictions, a) Comparison with QCD predictions calculated atNLO byFre- 
und et al. [24] based on MRST 2001 and CTEQ6 PDFs. b) Comparison with the colour dipole 
predictions of Donnachie and Dosch [28] and Favart and Machado with [31] and without [29] 
the DGLAP evolution of the saturating dipole (indicated as BGBK). The band associated with 
each prediction corresponds to the uncertainty on the measured t-slope. 
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